Turk J Phys 

27 (2003) ,1 - 10. 

© TUBiTAK 



Decoherence in two coupled qubits 

Kristian RABENSTEIN and Dmitri V. AVERIN 

Department of Physics and Astronomy, Stony Brook University, SUNY 
Stony Brook, NY 11794, U.S.A. 



Received February 2, 2008 



Abstract 

We have developed quantitative description of quantum coherent oscillations in the system of two 
coupled qubits in the presence of weak decoherence that in general can be correlated between the two 
qubits. It is shown that in the experimentally realized scheme of excitation of the oscillations, their 
waveform is not very sensitive to the magnitude of decoherence correlations. Modification of this scheme 
into potentially useful probe of the degree of decoherence correlations at the two qubits is suggested. 
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1. Introduction 

Despite the large number of successful demonstrations of quantum coherent oscillations in individual 
[1, 2, 3, 4, 5, 6] and coupled [7, 8] Josephson-j unction qubits, quantitative understanding of the details of 
these oscillations is so far quite limited. One of the most important open problems is decoherence, many 
aspects of which still remain to be understood. The purpose of this work is to develop theoretical description 
of decoherence in dynamics of a system of two coupled qubits. The approach we are using is based on the 
evolution equation for the density matrix in the Markovian approximation that is standard for description of 
weak decoherence. Although there are indications from the single-qubit experiments that more sophisticated 
approaches are needed for quantitatively accurate description of decoherence, the result we obtain within 
the simple scheme can be useful as the benchmark for more elaborate models. 

Motivation for studying decoherence in coupled qubits is provided by the recent experiment [7] with two 
coupled charge qubits, where it was found that the decoherence rate for quantum coherent oscillations in 
two qubits at the optimal bias point is with good accuracy factor-of-4 larger than the decoherence rate in 
effectively decoupled qubits. An interesting question for theory is whether this factor-of-4 increase of the 
decoherence rate is a numerical coincidence, or it reflects some basic property of the decoherence mechanisms 
in charge qubits. As will become clear from the discussion below, the theory developed in this work favors 
"numerical coincidence" point of view. Other aspects of decoherence in coupled qubits has been studied 
before numerically in [9, 10, 11]. 

In general, it is well understood that decoherence rates of different states of two coupled qubits can 
be quite different if the random forces created by the qubit environments responsible for decoherence are 
completely or partially correlated at the two qubits. Most importantly, in the case of complete correlation, 
the qubit system should have a "decoherence-free subspace" (DFS) spanned by the states |01), |10) [12, 
13, 14], since completely correlated external environments can not distinguish these states. In contrast, 
the subspace spanned by |00) and experiences decoherence that is made stronger by the correlations 
between environmental forces acting on the two qubits. So the role of the quantitative theory in description of 
decoherence in the dynamics of coupled qubits is to see to what extent subspaces with different decoherence 
rates participate in the qubit oscillations for different methods of their excitation. 
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2. The model and environmental correlations 

The Hamiltonian of the system of two qubits coupled directly by interaction between the basis-forming 
degrees of freedom (i.e., electrostatic interaction through finite coupling capacitance for charge qubits, or 
magnetic interaction for flux qubits) is: 

i=i,2 

where ct's denote the Pauli matrices, is the qubit interaction energy, Aj is the tunnel amplitude and ej 
is the bias of the jth qubit. Four energy levels of the Hamiltonian (1) are shown schematically in Fig. 1 as 
functions of the common bias e = ei = £2 of the two qubits. In this work, we will consider quantum coherent 
oscillations in the qubits biased at the "co-resonance" point [7], where £1 = £2 = 0. Such bias conditions are 
optimal for the oscillations. 

It can be shown explicitly that the occupation probabilities of the qubit basis states (that are of interest 
for us) are insensitive to the phases of the qiibit tunnel amplitudes Aj , so without the loss of generality we 
will assume that A^ 's are real. The Hamiltonian (1) at the co-resonance reduces then to 

Ho=Y.Aj4^^+^a^^->a(^K (2) 

i=i,2 

In the basis composed of eigenstates of the ai^^ operators, the Hamiltonian (2) can be diagonalized easily. 
Eigenenergies and eigenstates are: 

El = Q, 1^1) = ^[(7 + V){m + 111)) + (7 - r?)(|01) + |10))], 

E2 = -n, |V-2) = i[(r7-7)(|00) + |ll)) + (7 + r;)(|01) + |10))], (3) 
E3 = e, 1^3) = l[i^ + /3)(|00) - 111)) + (a - /3)(|10) - |01))], 
E, = -e, 1^4) = - «)(|00) - 111)) + (a + /J)(|10) - |01))], 

where 

and A = Ai + A2, 5 = Ai — A2. The states \kl) with {fc, 1} — {0, 1} in Eqs. (3) are the eigenstates of the 
operators cri^'^'' in the natural notations: (j^z''\kl) = {—l)''\kl) and G^z''\kl) = (— l)'|fc^). 

We assume that external environments responsible for the decoherence couple to the basis-forming degrees 
of freedom of the qubits. The interaction Hamiltonian is then: 

Hi=Y.m<^'^JK (4) 

The random forces S,i,2{t) acting on the qubits are in general correlated. To describe the weakly dissipative 
dynamics of the system in the basis of states (3) induced by the interaction (4) with the reservoirs, we use 
the standard equation for the evolution of the qubit density matrix p in the interaction representation (see, 
e.g., [15]): 

p = -[ dT{[m{t),[Hi{T),p]]), (5) 

J —00 

where angled brackets denote averaging over the reservoirs. Proceeding in the standard way, we keep in Eqn. 
(5) only terms that do not oscillate as functions of time with large frequencies on the order of eigenenergies 
of the system, and therefore lead to changes in p that accumulate over time. Equations for the matrix 
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0.0 

common bias 8 

Figure 1. Schematic structure of the energy levels of the two coupled qubits as functions of the common bias of 
the qubits. The inset shows the diagram of the decoherence-induced transitions between the levels at "co-resonance" 
point where the bias vanishes. 



elements pnm, n, m = 1, ...,4, of p in the basis (3) are transformed then as follows: 



fe 

+ E PkK^nk^'liliPn'i^l - ^m) + F;,^{ei - e„))] . 
(fe,0 



(6) 



Here Unm denote the matrix elements {n\a^J^\m), the last sum is taken over the pairs {k,l) of states that 
satisfy the "resonance" condition: 



Cfe - = En - Efi 



{k, I) ^ (n, m) , 



and 



The first term in Eqn. (6) represents "pure dcphasing" that exists when the system operators that couple 
it to the environment have non- vanishing average values in the eigenstates. As one can see explicitly from 
Eqs. (3), the average values ann of a^z^ are vanishing in all states, so that there is no pure dephasing term 
in the evolution of the density matrix in our case. The fact that ann are vanishing can be related to the 
vanishing slope of the system energies with respect to variations of the bias in the vicinity of the co-resonance 
point - see Fig. 1. Since all coefScients in the eigenfunctions (3) are real, the matrix elements Gnm are also 
real. For real Unm, imaginary parts of the noise correlators Fjij in the second term on the right-hand-side 
of Eqn. (6) represent the decoherence-induced shifts of the system energy levels. These shifts do not affect 
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decoherence and we will neglect them in our discussion. With these simplifications, Eqn. (6) takes the form 

Pun. = J2 [-iPnn^mJ2('^l^y£^ReF,,,{en-ek)+all^aiiReF,J,{ern-e,)) 

k 

k {k,l) 

where 

/oo 
dt(e,(i)e,-,(0))e-*. (8) 
-oo 

The function F12 characterizes correlations between the environmental forces acting on the two qubits. 
For instance, if the two qubits interact with different environments and ^1, ^2 arc uncorrolatcd, F12 = 0, 
whereas F12 = -Fii = F22, if the qubits are acted upon by the force produced by one environment coupled 
equally to the two qubits. While the correlators Fn and F22 are necessarily real, F12 can be imaginary, and 
F21 = Fi2. Non- vanishing imaginary part of F12 corresponds to the non- vanishing commutator [^1,^2] and 
implies that the two qubits are coupled to the two non-commuting variables of the same reservoir. While 
this is probably not very likely for qubits with the basis-forming degrees of freedom of the same nature 
(which in a typical situation should be coupled to the same set of environmental degrees of freedom), the 
non- vanishing ImFi2 should be typical if the qubits have different basis-forming variables. Using the spectral 
decomposition of the correlators Fjjt {u) and Swartz inequality, one can prove (similarly to what is done in a 
different context of linear quantum measurements [16]) that for arbitrary stationary reservoirs the correlators 
satisfy the inequality that imposes the constraint on ^12(0;): 

i=^ii(w)F22(a;) > |Fi2(a;)p. (9) 

If the reservoirs are in equilibrium at temperature T, the correlators satisfy also the standard detailed balance 
relations: 

Fjj,{-io)^e-^Fj,j{^). (10) 

Equation (7) with the noise correlators (8) govern weakly dissipative time evolution of the two coupled 
qubits in a generic situation. Below we use them to determine decoherence properties of quantum coherent 
oscillations of the qubits. Before doing this, however, we would like to briefly discuss applicability of our 
approach to realistic Josephson-junction qubits. As we saw above, one of the main features of Eqn. (7) is 
that the pure dephasing terms disappear at the co-resonance point and remaining decoherence is related to 
transitions between the energy eigenstates. This implies that within the approach based on Eqn. (7), the 
decoherence rates are on the order of half of the transition rates, whereas experiments with charge qubits (see, 
e.g., [4]) indicate that decoherence rates are larger than the transition rates even at the optimum bias point 
when the pure-dephasing terms should disappear. Apparently, this is related to the low-frequency charge 
noise [4, 17] that is coupled to qubit strongly enough for the lowest-order perturbation theory in coupling 
(5) to be insufficient. This implies that the theory presented in this work might be only qualitatively correct 
for realistic charge qubits, for which one should develop more accurate non-perturbative description of the 
low-frequency noise to achieve quantitative agreement with experiments. Our simple perturbative approach, 
however, should be applicable to flux qubits, where the low-frequency noise should not be as strong as in 
the charge qubits. 



3. Quantum coherent oscillations in coupled qubits 

One of the most direct ways of excitation of quantum coherent oscillations in individual or coupled qubits 
that will be discussed in this work is based on the abrupt variation of the bias conditions [1, 7]. If the qubits 
are initially localized in one of their basis states, e.g. )00), and abrupt variation of the bias brings them 
to the co-resonance point, the probabilities for the qubits to be in other basis states start oscillating with 
time. In the simplest detection scheme (realized, for instance, in experiment [7]) the probability for each 
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qubit to be in the state |1) is measured independently of the state of the other qubit. Quantitatively, these 
probabilities pj are obtained from the projection operators Pf. 

Finding explicitly the matrix elements of Pj from the wavefunctions (3), one gets: 

Pi{t) = 1 + {<^V + /37)Re[e-'--*(p42 W - PisW)] + {aj - /3r?)Re[e-*-+*(pi4W + P32(0)] , (H) 

P2{t) = ^ - (a7 + /3r?)Re[e-^'^-*(pi3(t) + P42(i))] + ("^ - 0'y)Re[e-'^+' {pi4{t) - P32 W)] , (12) 

where oj± = fi ± e, and as in the Eqn. (7), the matrix elements of the density matrix are taken in the 

interaction representation. Equations (12) and (7) show that the waveform of the coherent oscillations in 
coupled qubits is determined by the time evolution of the two pairs of the matrix elements of p: 

Pl3 = -Eispis + U-P42 , P42 = -r42P42 + W+PlS , (13) 
Pl4 = -ri4Pl4 + «-P32 . P32 = ~r32P32 + W+P14 • (14) 

The decoherence rates in these equations are determined by the rates of transitions between different 
energy eigenstates: 

ri3 = i(rl+^ + r(-) + r(+) + ri+)) , r,4 = i(r(-) + r(+) + r(+) + 4+^) , 

r32 = I (r(-) + r(-) + ri") + ri+)) , r42 = ^ (r(-) + ri") + r(+) + ri")) . (15) 

where labeling of the transitions is indicated in the inset in Fig. 1. Transition rates are: 

id' 

r(±) = Re^F.y (±a;_)4iVi{) , Pf ' = RcJ^ F,,,{±u;^)a^J^ ag'> . (16) 

The superscripts ± refer here to transitions in the direction of decreasing (+) or increasing (-) energy. Finding 
matrix elements anm from the wavefunctions (3) we see explicitly that transitions between the states 1 and 
2, as well as 3 and 4 are suppressed, since the corresponding matrix elements are zero, and that the rates 
(16) are: 

5 A — ^ ^ ,1^ 
^^) + Refr2(^--), 

^ A — ;y2 ^ ^ „ ,u 
_^) + Rei.,,(- + -), 

^^)-ReF,2(^ + -), (17) 
^^-^VReF,2(^-^). 



ri = 




(5A + 1/2 


^F22(l + 


r2 = 




5A + ;y2 


^^22(1 - 


r3 = 


^Fn(l + 


6 A + ;y2 

en ) + 


^F22(l- 


r4 = 


^Fn(l- 


(5A + 1/2 

en 


^i^22(l + 



en ' "n e 

The transfer "rates" w, v in Eqs. (13) and (14) are: 



Explicitly: 
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Equations (17) and (19) do not show the frequency dependence of noise correlators, which is the same, 
respectively, as in the Eqs. (16) and (18). 

Each pair, (13) and (14), of coupled equations can be solved directly by diagonalization of the matrix 
of the evolution coefficients with a non-orthogonal transformation. In this way we obtain for the pair of 
equations (13): 

Pi3(i) = ^—^ [pi3(0)(n+«_e-^+* + c^e-^-*) + c«_/942(0)(e-T+* - e"^"*)] , 

Pi2{t) = ^— ^ [p42(0)(u+«_e-^-* + c2e-^+*) + cw+pi3(0)(e-^+* - e"^"*)] . (20) 

where 

7± = (ri3 + r42)/2 ± [(ri3 - r42) V4 + u+u-] ^'^ , c = (ri3 - r42)/2 - [(Eis - E42) V4 + u+u-] ^'^ , (21) 

and /9i3(0), /942(0) arc the initial values of the density matrix elements that depend on preparation of the 
initial state. If, as in the experiment [7], the qubits are abruptly driven to co-resonance maintaining the 
state 1 00), these initial values are: 

Pi3(0) = J(7 + ni){a + 13) , P32(0) = J(7 - V){a + (3) , 

Pi4(0) = J(7 + '?)(/3-a), P42(0) = J(7 -»?)(«-/?). (22) 

Another type of initial conditions that will be discussed in this work is starting the oscillations from the 
state 1 10). In this case: 

Pi3(0) = J(7 - V){a - P) , P32(0) = J(7 + V){a - P) , 

Pi4(0) = J(7 -»?)(" + /?), P42(0) = J(7 + »?)(« + /?). (23) 

Equations (22) and (23) follow directly from the wavefunctions (3): Pnm(O) = {n\i){i\m), where \i) is the 
initial state. 

Solution of the other pair (14) of coupled equation is given by the same Eqs. (20) and (21) with obvious 
substitutions: u± —^ v±, T13 T14, T42 ^32- In this work, we are mostly interested in the low- 
temperature regime T <^ e,Q, when transitions up in energy can be neglected. In this regime, U-,V- — > 0, 
and equations for the evolution of the density matrix elements are simpUfied. For instance, for U- — > 0, 
c ~ ■u+u_/(ri3 — r42), and Eqs. (20) are reduced to: 

Pi3(i) = Pi3(0)e-^-* , p42{t) = P42(0)e-r-* + Pi3(0)(e-^^='* - 6"^-*) , (24) 

i 13 — i 42 

where now ri3 = (Fi +T2 + T4)/2 and r42 = F3/2. 

Time evolution (24) of the density matrix elements together with the rates (17) and (19), and initial 
conditions (22) and (23) determines the shape of the coherent oscillations in two coupled qubits. In the next 
Section, we discuss this shape in several specific situations. 



4. Results and Conclusions 

The shape of coherent oscillations determined in the previous Section depends on the large number of 
parameters: temperature, degree of asymmetry of qubit tunnel energies and couplings to the environments, 
frequency dependence of the decoherence, and strength and nature of the correlations between the two 
reservoirs. We analyze some of these dependencies below. 
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4.1. Experimentally- motivated case 

We start by considering the situation that is close to the experimentally realized case of oscillations in 
coupled charge qubits [7] . As we argued above, the correlations between environments in this case should be 
real: lmFi2 = 0. The oscillations are excited by driving the system to co-resonance in the initial state |00). 
Equations of the previous Section give in this case the following expression for the shape of the oscillations: 

Pi{t) = l-l [^e-^-* + S(e-r-* - -J^(e-r-* - 6"^-*))] cosa;_t 

Z » i 13 — i 42 

_i [(7e-r32« + £)(e-ri4« + _.^±_(e-ri4t _ 6-^32*))] cosw+i , (25) 



14 — J- 32 



where 



^_^_^JA + z/2 p V ^_ •|^_(_^^ + '^^_|_^_(_^ C-\ ^^^^^ ^ ^'^ D — 1 ^^^^"^ ^ 
ef2 n e ' ef2 fi e ' efi f2 e ' ef2 fi e 

Equation for P2(0 is the same with signs in front of (5, m+ and reversed. As a simplifying assumption 
we take F\\ = F22 = F. The decoherence rates in Eqn. (25) then are: 

ri3 = Fi.,)il - ^) + iF(._)(l + ^) + 1f.(._)(^ + ^) , u, = -i^(.,)§ , 

ri4 = F(._)(l + ^) + i^(.,)(l - ^) + \f.mC^ - ^) , = -F(._)§ , (26) 

r32 = \Fi.,){i - ^) + \F.i.,r- - ^) , r42 = ii^(.-)(i + ^) - i^e(..-)(^ + ^) , 

where Fc{u)) = ReFi2(a;). To enable comparison of these rates to those of individual qubits, we note that 
the rate of decoherence of oscillations in a single qubit with vanishing bias is equal to F(Aj)/2 for the jth 
qubit. 

The functions Pi,2(i) for the qubit parameters, 5, v, and F, close to those in experiment [7] are plotted in 
Fig. 2 under additional assumption that the decoherence is the same at two frequencies, uj^ and (The 
decoherence strength F was taken from the data for the single-qubit regime in [7].) The curves are plotted 
for the two situations, when decoherence is completely uncorrelated (F12 = 0) and completely correlated 
{F12 = F) between the two qiibits. One can sec that the difference between the two regimes is not very big 
numerically, with correlations between the two reservoirs leading to the effective decoherence rate that is 
increased in comparison with the uncorrelated regime by roughly 30 50%, although the description with a 
single decoherence rate is not quite appropriate quantitatively - sec Eqs. (25) and (26). 

The increase of the effective decoherence rate by correlations illustrated in Fig. 2 can be related to the fact 
that the initial qubit state, |00), belongs to the subspace where the correlations increase the decoherence rate, 
despite the mixing of this subspace with the DFS where the decoherence rate is decreased in the eigenstates 
(3) of the coupled qubit system. (Here and below we use the term "DFS" for the subspace spanned by the 
|01) and 1 10) states, although for interacting qubits it, strictly speaking, does not fully have the properties 
of real DFS.) This implies that increase of decoherence rate by correlations should be to a large extent 
insensitive to qubit parameters. This conclusion is supported by the case of identical qubits ((5 = 0), when 
u+ = v+ = and Eq. (25) is reduced to a very simple form: 

Mt) = l-\ [(1 + cosuj-t + (1 - ^)e-^-* cosu+t] , (27) 

and P2{t) = Pi{t). One can see from Eqs. (26) that both decoherence rates relevant for Eq. (27), ri3 and 
r32 increase with increasing correlation strength Fc- Equation (27) shows also that the description of the 
oscillation decay with a single decoherence rate can be quite inaccurate: for weak interaction, u < fl the 
amplitudes of the two (high- and low-frequency) components of the oscillations are nearly the same while 
their decoherence rates can be very different. 
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Figure 2. Probabilities pj to find jth qubit in the state |1) in the process of quantum coherent oscillations starting 
with the state |00) of two coupled qubits. The decoherence strength is F = 0.08A. Solid and dashed lines correspond, 
respectively, to the decoherence that is uncorrelated {Fc = 0) and completely correlated (Fc = F) between the two 
qubits. 



4.2. Excitation into the DFS 

Now we discuss decoherence properties of the oscillations in coupled qubits in the case when they start 
with the initial qubit state |10). We note that in the case of experiment similar to [7], such an initial 
condition would require separate gate control of the two qubits, since the bias change bringing them into 
CO- resonance is different in this state for the two qubits. Since the state |10) belongs to the DFS in the case 
of completely correlated noise, one can expect that oscillations with these initial conditions will be more 
sensitive to the degree of inter-qubit decoherence correlations than oscillations with |00) initial condition, 
and that the effective decoherence rate will decrease with correlation strength. All this indeed can be seen 
from Eqs. (24) with the initial conditions (23) that correspond to the |10) state. Under the same assumptions 
as were used in Eqn. (25), we get for the now different pi{t) and P2{t)- 

PAt) = \ - ^{ [^.e^^-* + B,(e-i^-* + izlZ^(e-r.3* _ e-^^-^)] cosu._i 
la i 13 — J- 42 

+ [C,e-^-* + ^.(,-r.t _ izi)!^(e-ri.* _ e'^-t))] ^ = 1^2, (28) 

i 14 — i 32 



where 



>li = 1 + TT- + + - > Si = 1 ■ 



e 



and the amplitudes A2, B2 ,€'2 , -D2 are given by the same expressions with d ^ —S. 
For identical qubits Eqn. (28) reduces to: 

Pjit) = 1- [(1 + cosw.t + (1 - ^)e"'^"* cosa;+i] . (29) 
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This expression and Eqs. (26) show that in contrast to Eqn. (27), the decohercncc rate of the low- frequency 
component that has larger amplitude is strongly suppressed by the non-vanishing inter-qubit noise corre- 
lations Fei T42 = 1(1 + i'/n)[F{uj-) - Fc(w_)]. This means that the shape of the coherent oscillations in 
coupled qubit starting with the state |10) should indeed be more sensitive to the strength of these corre- 
lations than the shape of the oscillations starting with the |00) state. As one can see from Fig. 3 which 
shows the shape (28) of the |10) oscillations for the same set of experimentally realized parameters as in 
Fig. 3; this conclusion also remains valid in the case of not fully symmetric qubits. Even in this case there 
is a pronounced weakly decaying component of the oscillations if the decoherence is completely correlated 
between the two qubits. For partial correlations, the effective decoherence rate is reduced. 

0.8 



d! 0.4 



0.0 

0.0 5.0 10.0 

1.0 



d: 0.5 



0.0 

0.0 5.0 10.0 

Figure 3. Probabilities pj to find jth qubit in the state |1) in the process of quantum coherent oscillations starting 
with the state |10) of two coupled qubits. Qubit parameters are the same as in Fig. 2. Solid, dotted, and dashed lines 
correspond, respectively, to the decoherence that is uncorrelated {Fc = 0), partially {Fc = 0.5F), and completely 
{Fc = F) correlated between the two qubits. 

In summary, we have developed quantitative description of weakly dissipative dynamics of two coupled 
qubits based on the standard Markovian evolution equation for the density matrix. This description shows 
that decoherence properties of currently realized oscillations in coupled qubits are not very sensitive to 
inter-qubit correlations of decoherence, while relatively simple modification of the excitation scheme for the 
oscillations should make them sensitive to these correlations. 
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